using LinearAlgebra
using MKL
BLAS.set_num_threads(1); println(BLAS.get_config())
const function_path = pwd(); #"/path/to/dir/"
const save_dir = pwd(); #"/path/to/dir/"
const potential_dir = pwd(); #"/path/to/dir/"
const pstring = "depth2.446-nlvl5-phi_echo1.5707963267948966-beta4.0-dr50.-rmax400.-astig106.4-dchi0.-gam_motion60.0"
const num_lvls = 5
const max_ho_level = 40
const aberration_type = "-dr50.-rmax400.-astig0-dchi0."
const Omega_pulse = 38000.0
const Omega_drive = 5180.0
const phi_drive = 4.71238898038469
const phi_echo = 1.5707963267948966
const beta = 4.0
const trap_depth = 2.446
const gam_deph = 12.5
const gam_deph_motion = 60.0
const krylov_tol = 1.0e-10
include("$function_path/nacs-tweezer-dvr.jl")
include("$function_path/nacs-compute-interactions.jl")
include("$function_path/onebody-drive-dephasing.jl")
include("$function_path/onebody-echo-dephasing.jl")
include("$function_path/onebody-xy8-dephasing.jl")
